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Abstract 

In this paper we report comprehensive experimental and theoretical investigation of magnetic 
and electronic properties of the intermetallic compounds Pr2Fei7 and Gd2Fei7. For the first time 
electronic structure of these two systems was probed by optical measurements in the spectral 
range of 0.22-15 //m. On top of that charge carriers parameters (plasma frequency £1 and relax- 
ation frequency 7) and optical conductivity <j(ui) were determined. Self-consistent spin-resolved 
bandstructure calculations within the conventional LSDA+U method were performed. Theoretical 
interpetation of the experimental cr(u) dispersions indicates transitions between 3d and 4p states of 
Fe ions to be the biggest ones. Qualitatively the line shape of the theoretical optical conductivity 
coincides well with our experimental data. Calculated by LSDA+U method magnetic moments per 
formula unit are found to be in good agreement with observed experimental values of saturation 
magnetization. 

PACS numbers: 71.20.-b, 78.20-e, 75.50.Ww 
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I. INTRODUCTION 



During last two decades Fe-rich intermetallic systems stays under investigation from 
both experimental and theoretical points of view because its anomalous magnetic proper- 
ties. At present the intermetallic compounds with rare-earth elements draw a significant 
interest because of their well pronounced Invar properties and also as possible new cheap 
high-energy storing materials for permanent magnets . 1 ! 2 ! 3 ! 4 ! 5 ! 6 ' 7 Rare-earth intermetallic com- 
pounds R,2Fei7 (R - a rare-earth ion) are well-known due to large magnetic moment of Fe 
ions and yet their comparatively low Curie temperature Tq- Further increase of the Tc 
of the compounds can be achieved by implantation of interstitial atoms^ 2 or introducing 
a low content of non-magnetic impurities substituting Fe atoms . 3 ! 4 ' 5 ! 6 ' 7 In such extended 
systems, compared to the parent R2Fei7 compounds, upon the doping Tc grows more than 
twice in some cases. This increment of Tq was explained primarily as a result of the lat- 
tice expansion with substitution of Fe ions for larger radii ions. According to this concept, 
Fe-Fe interactions are ferromagnetic or antiferromagnetic for interatomic distance larger or 
smaller then the critical value 0.25 nm.— However, a number of experiments^ ' 10 ' 11 showed 
that in the substitutional R2Fei7_ x Si x alloys the crystal lattice contracts and the magnetic 
moment per unit volume decreases, while Tc grows. Later experimental facts are in appar- 
ent contradiction with the model.— Clearly, a simple approach based on distance-dependent 
exchange interaction is not sufficient to explain the changes of magnetic properties of these 
compounds with substitution of Fe for non-magnetic ions. 

Past life several results of band structure calculations of the systems from the family 
are available. For instance the spin-polarized calculations of electron spectra of Y 2 Fei7N 3 ,— 
SmaFe^-xMz (M=A1, Ga, Si)^ and Nd 2 Fei 7 -xM x (M = Si, Ga) 14 compounds showed that 
non-magnetic impurities modify the shape and width of the spin-up and spin-down densities 
of states N(E). Jaswal et ah— suggested to explain the increase of T c in Y 2 Fei 7 N 3 via 
changes of the density of states at the Fermi level N(Ep) owing to the lattice expansion 
by means of the Mohn-Wohlfarth static spin-fluctuation model.— An analysis of the ex- 
perimental optical , 7 ' 16 ' 17 low-temperature heat capacity 7 and photoemission 1 ^ data of some 
pseudobinary alloys K2Fen- x M x with M = Al, Si revealed a qualitative correlation between 
the Curie temperature and parameters of electronic structure in frame of such an approach. 
Recent experimentally measured optical conductivity of Ce2Fei7 was interpreted in terms of 
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band structure obtained within local density approximation (LDA).— 

In accordance with the foregoing systematic study of the electronic structure of R2Fei7 
compounds and their modifications with substituted Fe ions further investigations are of 
fundamental importance. In this work theoretical calculations of the electronic properties 
together with experimental magnetic and optical measurements were performed. The paper 
is organized as follows: in the section |H] experimental details and results for Pr2Fei7 and 
Gd2Fei7 are presented. For instance, subsection III Al is devoted to sample preparation, X-ray 
diffraction structural analysis and magnetic measurements conditions. Subsection III Bl pro- 
vides description of optical experiments. Section ITTTl contains results of LSDA+U 20 computa- 
tions of electronic structure and magnetic properties of Pr-Fe and Gd-Fe systems. Structure 
of experimentally observed optical conductivity curves is anatomized in the section II VI At 
the end we briefly summarize our paper with the section 

II. EXPERIMENT 

A. Samples and magnetic measurements 

The compounds Pr 2 Fei7 an d Gd2Fei7 were prepared by induction melting in alumina 
crucible under argon atmosphere. The ingots were homogenized in the high-purity argon 
atmosphere at ~1300 K. The purity of the alloys was checked using standard X-ray diffrac- 
tometry in Cu Ka radiation. The samples were found to be polycrystalline, single-phase 
and have rhombohedral structure of the T/^Zn^-type^i (space group R3m) for Pr 2 Fei7 an d 
hexagonal crystal structure of the Th2Nin-type2& (space group P6s/mmc) for Gd2Fei7. The 
measured lattice parameters a and c are close to those published earlie r 21 i 22 and are shown 
in the Table UJ Spherical specimens of 2-3 mm in diameter were used for magnetic measure- 
ments. For following optical studies the specular surface of the samples was prepared by 
mechanical polishing with diamond pastes. 

The Curie points of the compounds were determined from temperature dependencies 
of ac susceptibility. The ac susceptibility was measured by a differential method in an ac 
magnetic field of 8 Oe with a frequency of 80 Hz. The saturation magnetizations at T=4.2 K 
were determined from the isothermal magnetization measurements carried out by means of 
vibrating sample magnetometer in magnetic fields up to 20 kOe. The magnetic parameters 
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TABLE I: Experimental structural, magnetic, and electronic parameters of the intermetallic com- 
pounds Pr2Fei7 and Gd2Fei7: lattice parameters a, c; unit-cell volume at room temperature V; 
Tc - Curie temperature; M s - spontaneous magnetization at T = 4.2 K; relaxation 7 and plasma 
frequencies; N e ff - effective concentration of conduction electrons. 



Compound 


a (A) 


c (A) V (A 3 ) T c (K) M s (/is/f-u.) 


7 (lO^s" 1 ) 


n 2 (io 30 s- 2 ) N eff 


(10 22 cm- 3 ; 


Pr 2 Fei 7 


8.579 


12.472 795.0 294 36.1 


1.9 


21.3 


0.67 


Gd 2 Fei7 


8.496 


8.341 521.4 466 21.2 


1.5 


19.5 


0.61 



of the compounds are listed in Table |U They are a little different from previously reported 
data . 23 i 24 The differences among authors may originate from some deviations of their samples 
from stoichiometry. 

B. Optical measurements 

Investigation of optical properties of Pr 2 Fei7 an d Gd 2 Fei7 were carried out at room tem- 
perature by ellipsometric Beattie technique.— Spectroscopic ellipsometry is based on the 
fact that the state of polarization of incident light is changed on reflection. The optical 
constants - refractive index n and absorption coefficient k - were measured in the spec- 
tral range of fta>=0.077-5.6 eV (u is a cyclic frequency of light) with accuracy of 2-4%. 
From n and k, the real ex(ou) = n 2 — k 2 and imaginary 62(00) = Ink parts of the com- 
plex dielectric function e(u), the optical conductivity <t(uj) = nku/2iT, and the reflectance 
R(u) = [{n — l) 2 + k 2 ]/[(n + l) 2 + k 2 } were derived. Measurements of reflection spectra 
followed by the Kramers-Kronig analysis were applied to determine the optical parameters 
in the short-wave range (Hlu— 5.6-8.5 eV). 

The results of optical study (n, k, ei, e 2 , R as functions of a;) for the Pr and Gd compounds 
are shown in Fig. It is seen that the optical properties of both compounds are rather 
similar. All directly measured dispersions and further derived quantities are characterized 
by the broad feature with maximum at the photon energies near 1 eV. As it follows from 
the lineshape of the e 2 (o;), there is a strong absorption region at energies ftw >1 eV. With 
increase of the wavelength to infrared range (A >2 jum) (inset of Fig. the nonmonotonic 
behaviour of n and k is changed by a smooth growth related to the domination of free-electron 
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absorption. Such a character of the frequency dependencies for the optical constants together 
with negative quantities e%, as a rule, is typical for the metal- like solids. As one can see 
from dispersion curves R{uj), in the low-energy range the reflectance exhibits the high values. 
The analysis of the energy dependence of e% and e 2 i n this spectral interval, corresponding to 
intraband electronic excitations, makes it possible to determine the plasma frequency Q and 
the effective relaxation frequency 7 of free charge carriers. Within the assumption that light 
absorption for these energies has Drude character, the parameters Q 2 = uJ 2 (e\ + e|)/ei and 
7 = e 2 uj/ei were computed. In the long- wave region A >8 /im, 7 and fl 2 become frequency 
independent. The values of 7 and Q 2 in this energy interval were used then to estimate the 
effective concentration of conduction electrons N e ff = f2 2 m/47re 2 (m and e are the mass and 
the charge of a free electron respectively). All these parameters obtained from experimental 
data treatment are presented in Table |U 

Optical conductivities cr(u) for both compounds are displayed in Fig. |21 A monotonic 
increment of the experimental cr(u) dispersion observed in the low-energy range is related 
to the Drude mechanism of electron excitations. For both systems intraband (Drude-like) 
contribution to the optical conductivity were computed according to the relation <Ji n tra{^) = 
f2 2 7 /47t(cj 2 + 7 2 ) and drown in Fig. [21 black dotted lines. Corresponding values of Q 2 and 
7 are given in Table |1J The magnitude of these contributions falls down sharply with 
energy and becomes insignificant above 0.5 eV. The values of static conductivity cr/ ntra (0) 
are estimated to be 0.89*10 16 s _1 for Pr2Fei7 and 1.03*10 16 s _1 for Gd2Fei7 correspondingly. 
With increase of photon energy the <j(uj) curves show a very intense asymmetric stuctures in 
the near infrared region of spectra at ~1.2 eV. These structures have a pronounced shoulder 
on high-energy side and abrupt low-energy edge. Such a behaviour is a typical manifestation 
of the predominance of the interband absorption in this energy interval. The similar shape 
of the a{uj) curves was early observed in Y 2 Fei 7 , Ce 2 Fei 7 , and Lu 2 Fei 7 compounds. 17 



III. ELECTRONIC STRUCTURE CALCULATIONS 

To calculate electronic structure of the intermetallic compounds under investigation 
LSDA+U method 20 within the TB-LMTO-ASA package (Tight Binding, Linear Muffin-Tin 
Orbitals, Atomic Sphere Approximation)^ was applied. Experimentally obtained values of 
lattice constants for both Pr and Gd systems given in Table were used in our calcula- 
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tions. Atomic spheres radii were chosen as R(Pt) = 3.91 a.u. and R(Fe) = 2.62 a.u. for 
Pr 2 Fei7 and i?(Gd) = 3.72 a.u. and R(Fe) = 2.66 a.u. for Gd 2 Fei7. Orbital basis consists 
of 6s, 6p, 5d, and 4/ muffin-tin orbitals for Pr or Gd and 4s, 4p, and 3c? for Fe sites. The 
calculations were performed with 32 irreducible k-points (6x6x6 spacing) in the first Bril- 
louin zone. Parameters of direct U and exchange J Coulomb interactions for Gd and Pr 
ions were calculated by constrained LDA method.— For Gd 2 Fei7 we obtained £/c?d=6.7 eV 
and Jc,i=0.7 eV (similar values were determined previously for elemental Gd2&), and for 
Pr 2 Fei7 Up r =4.9 eV and Jp r =0.6 eV. Corresponding values of U and J were applied to 
Gd and Pr compounds in frame of the LSDA+U method. To note, in the present work 
we do not take into account local Coulomb interaction on Fe ions since constrained LDA 
gives surprisingly large U values.— However, for elemental Fe it was shown that account of 
Coulomb correlations is important to describe semiquantitatively temperature dependence 
of magnetization.— 9 But on the other hand such approach reproduces main structures of 
LSDA DOS with slight modifications of the Fe 3d DOS and thus in our case will not affect 
resulting dispersions of optical conductivity strongly. 

Total magnetic moments per formula units obtained within LSDA+U method (33.79 hb 
for Pr 2 Fei7 an d 21.89 [ib f° r Gd 2 Fei7) are in good agreement with measured experimental 
data (see Table |1J). Values of local magnetic moments for different sites for Pr 2 Fei7 and 
Gd 2 Fei7 compounds are listed in Table |HJ Pr and Gd 4/-shells are computed to be fully 
polarised with local moments close to its ionic values while Fe ions have local moments values 
almost equal to its elemental Fe magnitude. It is remarkable that local moments on rare- 
earth sites during selfconsistent loops become oppositely directed to those on Fe sites. One 
should also mention that initial value of local moments on Pr and Gd ions was taken to be 
zero. Thus R and Fe sublattices in our calculations are obtained to be antiferromagnetically 
ordered. 

Calculated in the present work within the LSDA+U method partial Fe 4p, 3d and rare- 
earth 5c?, 4/ DOS for spin-up (f) and spin-down (j) projections of local spin moments for 
Pr 2 Fei7 and Gd 2 Fei7 are shown in Figs. El and EJ It is seen that main spectral weight is 
located in the Ep ± 5 eV energy range around the Fermi level. Structures of these DOS are 
rather similar for both compounds and are qualitatively close to the DOS of ferromagnetic 
elemental iron in bcc structure.— Both these spin-polarized DOSs have a two-peak structure 
related to the 3d states of Fe ions. The Fermi level Ep set to zero lies near the minimum 
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TABLE II: Calculated within the LSDA+U method values of local magnetic moments for different 
sites in Pr2Fei7 and Gd2Fei7 compounds. 





Pr 2 Fei 7 




Gd 2 Fei 7 


Site 


M(fi B ) 


Site 


M(fi B ) 


Pr(6c) 


-2.08 


Gd(26) 


-7.13 


Fe(6c) 


2.21 


Gd(2d) 


-7.20 


Fe(9d) 


2.23 


Fe(4/) 


2.31 


Fe(18/) 


2.13 


Fe(6 5 ) 


2.10 


Fe(18/i) 


2.34 


Fe(12j) 


2.40 






Fe(12fc) 


1.87 



between these two peaks in the spin-down channel and on the upper edge in the spin-up 
channel. "Spin-up" states of Fe ions are almost completely occupied while "spin-down" are 
nearly half-filled. The narrow intensive peaks at 3 eV (j) and 2.7 eV (T) (Pr 2 Fei7) a ls° & t 
—7.3 eV (|) and 3.5 eV (f) (Gd 2 Fei 7 ) belong to 4/ states of rare-earth ions. The intensities 
of Fe 4s, Ap and rare-earth 6s, 6p, and 5d contributions to the DOS are considerably smaller. 



IV. ANALYSIS OF OPTICAL CONDUCTIVITY STRUCTURE 

Calculated LSDA+U band structures presented in the preceding section were used to 
interpret experimental optical conductivity <j(oj) for both intermetallic systems under con- 
sideration. In order to calculate theoretical optical conductivity u t h eor {uj) and anatomize its 
ingredients we applied rather simplified technique similar in spirit to Ref.^. Namely, we 
computed following convolutions representing all possible optical transitions 

= -jrj J E N u( E ) N u'( E - (i) 

where N^(E) and N^,(E) are partial DOSs of the same ion i with the same spin s and 
orbital quantum numbers I and /' shown in Figs. H3 and 0] The orbital quantum numbers are 
related by dipol selection rule I — I' = ±1. Total theoretical optical conductivity is linear 
combination of different contributions (Q): 

Otheor{u) = ^ miS ^'( U; )' ( 2 ) 
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where rrii is the degrees of site degeneracy (Wyckoff positions). 

Figs. 121 and El display corresponding a t heor{u) as solid lines. The overall shape of optical 
conductivity dispersion curves for both compounds exhibits a broad structure with the 
peak at ~2 eV. On the whole the theoretical calculations reproduce main features of the 
experimental ct{uj) curve (Fig. |2J): (i) the width of the intense absorption band with distinct 
maximum, (ii) the sharp threshold in the 0.5-1 eV range and (iii) the gradual diminution 
of the high-energy slope. 

At the same time there are certain discrepancies. The theoretical peaks are slightly shifted 
towards the higher energies in contrast to the experimental ones. Also high-energy contri- 
butions at 4 eV are slightly overestimated. It may have following reasons. As was shown by 
Lichtenstein et alM inclusion of local dynamical Coulomb interactions into consideration in 
elemental Fe will lead to about 30% of correlation bandwidth narrowing. In our case DOS in 
the vicinity of the Fermi level consists mostly from Fe 3c? states which are almost identical to 
those of elemental bcc iron. Roughly one can say that for our intermetallic compounds it can 
slightly improve comparison with experiment. First of all peak at 2 eV could be moved on 
0.5-0.6 eV towards lower energies because of correlation narrowing. Second spectral weight 
at 4 eV can be transfered to lower energies and it can give high-energy tail of theoretical 
optical conductivity closer to experimental one. Furthermore for both systems <j t heor{u) 
shows a significant interband absorption in the low-energy range (hw <0.6 eV) which was 
not confirmed by the measurements. It is possible that the low-energy absorption, predicted 
in theory, may be partially disguised in the experimental <j(oj) curves by the strong Drude 
rising. The Drude contribution was not taken into consideration in our theoretical model. 

Densities of states near the Fermi energy for both compounds mostly consist of the 3d- 
states of various crystallographically inequivalent iron ions (see Figs. H3 and 0J). A line shape 
of our LSDA calculated Fe ions DOS is very similar to the known DOS of elemental Fe 
in the bcc {hep) phase .'^ 2 For the bcc iron it was shown^ that the optical conductivity 
determined using the Berglung-Spicer approach^ agrees rather well with theoretical calcu- 
lations accounting matrix elements^ and also with experimental results.-^ This fact enables 
us to suppose that matrix elements formalism does not play an important role in case of 
intermetallic compounds R^Fe^. In this work we can report our theoretical curves to be in 
a reasonable qualitative agreement with experimental data. 

To analyze the line shape of experimental <j(uj) we provide detailed description of different 
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contributions to a t heor(^)- As it is seen in Fig. El for both systems biggest contribution to 
Ctheor(^) comes from 3d-4p interband optical transitions for Fe ions through the spin-down 
channel. That gives a high absorption peak at ~2 eV (see dotted curves in Fig. EJ). These 
particular transitions mostly govern low-energy range behaviour below 1 eV as well. Second 
largest contribution is 3d-4p transitions for Fe ions but in the spin-up channel. In Fig. |5] 
one can see it as a rather broad structure with maxima at ~4 eV (dot-dashed curves). The 
contributions from Fe 4s-4p (Fig. dashed- dot- dot lines) and Pr(Gd) 4f-5d (Fig. El thin 
solid line) transitions are substantially smaller; and magnitude of convolutions of rare earths 
5d-6p and 6s-6p states is negligible. 

V. SUMMARY 

Comprehensive experimental investigation of structural, magnetic and optical properties 
of the intermetallic compounds Pr2Fei7 and Gd2Fei7 was performed during this work. Re- 
fined structural parameters are found to be in good agreement with previous data. Measured 
magnetic properties: (i) Curie temperatures Tc=294K (466K) and (ii) saturation magneti- 
zations (at T=4.2K) 36.1 (21.2) /^/f.u. for Pr-Fe (Gd-Fe) systems respectively also agree 
well with available in the literature experimental data. We also report measured for the first 
time optical constants n and k observed at room temperature by ellipsometric Beattie tech- 
nique in the spectral range of 0.22-15 /mi. This experimental data allows us to determine 
charge carriers parameters (plasma Q and relaxation 7 frequencies) and optical conductivity 
a(uj). 

To model magnetic and optical properties of Pr2Fei7 and Gd2Fei7 we did self-consistent 
spin-resolved calculations within the LSDA+U method. Calculated by LSDA+U method 
magnetic moments per formula unit describe well observed experimental values. Further- 
more experimental optical conductivity c(u) was interpreted in terms of convolutions be- 
tween partial densities of states for the same ion applying dipole selection rule for orbital 
quantum number. Overall the line shape of experimental and theoretical optical conductiv- 
ity curves was found to be qualitatively very similar. By anatomizing different contributions 
to (Jtheor(w) it was understood that transitions between 3d and 4p states of Fe ions give the 
biggest contribution. The intense peaks in optical conductivity below 1 eV and around 2 eV 
are predominantly formed by the transitions in the "spin-down" channel while the second 
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large contribution from "spin-up" channel transitions provides broad feature with the max- 
imum at about 4 eV. The other contributions from rare-earth Af-5d and iron As-Ap optical 
transitions are almost negligible. 
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FIG. 1: (Colour online) Experimental optical constants n and k (panel a), dielectric functions ei and 
€2 (panel b), and reflectivity spectra R (panel c) for Pr2Fei7 (full circles and solid line) and Gd2Fei7 
(empty circles and doted line) compounds. Doted and solid lines above 5.6 eV represent values 
obtained by Kramers-Kronig transformation from reflection spectra. Inset presents dependences 
of the optical constants on wavelength for an expanded view in infrared region. 
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FIG. 2: (Colour online) Experimental (dark circles and dashed lines) and calculated (solid lines) dis- 
persion dependencies of optical conductivity for Pr2Fei7 (a) and Gd2Fei7 (b) compounds. Dashed 
lines represent ct(uj) values obtained by Kramers-Kronig transformation from reflection spectra. 
Doted lines correspond to intraband contributions estimated by Drude formula. Calculated values 
are in arbitrary units. 
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FIG. 3: (Colour online) Spin-resolved total (a) and partial 5d and 4/ DOS of rare earth ion (Pr) 
(b) and 3d and 4p densities of states of Fe (c) for Pr2Fei7 compound obtained from the LSDA+U 
calculation. DOSs were smoothed using adjacent averaging in the interval 0.1 eV. The Fermi level 
corresponds to zero. 
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FIG. 4: (Colour online) The same as Fig.EJbut for Gd 2 Fei 7 . 
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FIG. 5: (Colour online) Calculated total (thick solid lines) and partial contributions in optical con- 
ductivity for Pr2Fei7 (a) and Gd2Fei7 (b) compounds. Dashed, doted, and dash-doted lines show 
total, spin-up, and spin-down contributions from 3d-4p transitions in Fe ions respectively. Dash- 
dot-doted lines represent contribution from 4s-4p transitions in Fe, and thin solid lines correspond 
to 4/-5d transitions in rare-earth ions. 
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